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Abstract 

O ■ 

■ The periodic solutions of a type of nonlinear hyperbolic partial differential equations with a 

/"^v ' localized nonlinearity are investigated. For instance, these equations are known to describe several 

acoustical systems with fluid-structure interaction. It also encompasses particular types of delay 
differential equations. These systems undergo a bifurcation with the appearance of a small amplitude 
periodic regime. Assuming a certain regularity of the oscillating solution, several of its properties 
' around the bifurcation are given: bifurcation point, dependence of both the amplitude and period 

with respect to the bifurcation parameter, and law of decrease of the Fourier series components. All 
the properties of the standard Hopf bifurcation in the non-hyperbolic case are retrieved. In addition, 

1 r*| this study is based on a Fourier domain analysis and the harmonic balance method has been extended 
to the class of infinite dimensional problems hereby considered. Estimates on the errors made if the 
Fourier series is truncated are provided. 
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! 1 Introduction 
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1.1 Physical motivations 
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The starting point of the present work is due to physical motivations. The first objective is to justify 
a thirty-year-old conjecture assumed in order to obtaiminformation on the bifurcations encountered in 
woodwind musical instruments. It is now established |Tjf2| that woodwind instruments are nonlinear 
systems having a transition between a steady state and an oscillating regime, for a certain value of the 
pressure in the mouth of the musician. These instruments possess a fluid-structure interaction which 
implies the presence of partial differential equation (PDE) along with a nonlinearity. Analytical calcu- 
lations have been done for small oscillations around the bifurcation, giving valuable information on the 
oscillation with respect to the parameters of the instrument. But these calculations are possible only if a 
conjecture, first presented by Worman IjS], is assumed. Suppose the instrument emits a periodic sound of 
small amplitude, and the acoustic pressure po inside the mouthpiece is oscillating with period r. Let p a 
be the atmospheric pressure (or any other suitable pressure scale) and define the dimensionlcss quantity 
P — Po/Pa- Let P{u)q) be the Fourier component of p associated with the angular frequency u) q — 2qn jr. 
Worman's conjecture states: 

|P(w,)| < forallqeZ*, 

where, since it is of small amplitude, |P(u>i)| < 1. This key hypothesises not obvious and arguments 
to justify it have been given in the theoretical work of Grand et al IjljT However, the Fourier series 
has been truncated in the calculations of this latter work and the control of the remaining tail was 
not investigated. The first aim of the present study is to show that the weaker assumption of a twice 
continuously differentiable solution for the model of woodwind instruments implies Worman's conjecture. 
This former regularity assumption, allowing to control the tail of the Fourier series, is quite natural from 
a physical point of view. 

In addition to this first motivation, the hypotheses on the equations are weak enough to hold for 
a broader class of physical systems. The present approach is especially suitable for vibrating systems 
described by equations in the frequency domain with the formalism of transfer functi^pn^. This study 
proposes an extension of the harmonic balance method, widespread in engineering [5V [dj , and of the 
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underlying Galerkin procedure, to infinite dimensional systems governed by partial differential equations 
with a localized nonlinearity. This generalization holds for the case of small amplitude solutions which 
are twice continuously differentiable with respect to time. The present study justifies in a mathematical 
way the frequency domain expression, the truncation of the Fourier series, and the analysis around the 
bifurcation threshold. 

1.2 Mathematical framework 

Properties of the solutions of a nonlinear system around a bifurcation point can be obtained djy general 
techniques and theorems such as the reduction to the center manifold and normal forms , or 
the Hopf bifurcation theorem 1 )1 Ij . Nevertheless, in the case of PDEs involving the wave equation, 
which makes them hyperbolic, the hypotheses required for the use of these powerful theorems may 
not be satisfied. In such systems, the norm of the resolvent associated to the linearized operator in a 
neighbourhood of a bifurcation may not decrease when the spectral parameter goes to infinity. This 
is the case for woodwind musical instruments. The present study shows that although the hypotheses 
of the Hopf theorem are not satisfied, several of the results given by this theorem can be retrieved, 
provided the periodic solution is of small amplitude and twice continuously differentiable. The existence 
of such a solution is an hypothesis of the present work. Notice however that several studies show that 
nonlinear systems with PDEs possess such tt^ u j^ r ^solutions, for example the case of the wave equation 
with nonlinear terms and damping, see e.g. |T27T37~14 ■ 

After having settled the hypotheses satisfied by the systems under study, a corresp onde ncy between 
expr ession s of the model in the time and the frequency domain is establish^din lemma Ib.lT Then t heo- 
rems IHTlfe ves required conditions for an oscillation to appear and theorems fa. 1 1 and lb .11 with corollarv lb.^l 
give the dependence of both the Fourier harmonics and the frequency on the bifurcation parameter, near 
its birth point. As expected, it is similar to the finite dimensional case and what is given by the finite 
dimensional Hopf bifurcation theorem. In addition, this study gives also a valuable estimation of the law 
of decrease of the harmonics of the solution with respect to the multiple of the fundamental frequency in 
theorem lb.il This justifies the process of truncation used in the harmonic balance method and gives an 
error term on it. The last section comes back to the physical motivations and is devoted to an application 
in the field of musical acoustics. It introduces the equations modeling a reed woodwind instrument and 
shows how the mathematical results can be applied to this example. 
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2 The system 

Let (eTcM denote the time and Id" the space of configurations, for n e N. Let X = L 2 (T x X) 
be the separable Hilbert space of square integrable real- valued functions on the domain T x X. Only 
time periodic solutions are investigated and Tx I is supposed to be a compact domain of R n+1 . The 
nonlinearity takes place at a particular location xo S X in space and the infinite dimensional dynamical 
system under study is described formally by 

Ap = S(x - X )BG('J,p), (1) eqdep 

where p G X is the state of the system, the real number 7 is tp b-c called the bifurcation parameter, A 
is a linear partial differential operator. The right hand side of fl]) models the localized nonlinearity with 
Q being a nonlinear function of t J>^ d B a linear operator acting on Q and the symbol 8 being the Dirac 



distribution. As a consequence, $T|) must be understood in the distribution sense. In addition, since Tx X 
is compact, some linear boundary conditions are associated to |ff . fn a woodwind musical instrument, 



for instance a clarinet, the function p represents the air pressure inside the bore (resonator) and Ap is 
the wave equation (which may contain dissipative terms) in this domain of finite length. The position 
xq is the top of the instrument, where the tip of the reed is located. This flexible part oscillates when 
the musician plays and influence the air flow entering the instrument. The parameter 7 is proportional 
to the pressure inside the mouth. The function Q is the air flow through the clarinet, it depends in a 
nonlinear manner on the pressure difference between the player's mouth (7) and the resonator (p). Here, 
B is the time derivative operator. It is necessary to blow hard enough in order to obtain a sound, in other 
words the dynamical system encounters a bifurcation for some particular value of 7. It is important to 
understand how the oscillating frequency of the sound depends on the parameters of the instrument and 
the instrumentalist. Form example, obtaining this frequency will help understand how a beginner should 
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act to avoid a high frequency squealing oscillation. Since A Ls ^ differential operator one can always find a 
solution p, well defined as a distribution (see e.g. Schwartz ITS' section VII-10 and the examples therein). 
Let p be the Green function (distribution), solution of Ap = 8{x — £)8(t) and the boundary conditions. 
Thus p is given by the convolution p = p * 8{x — xq)BQ( , j ,p). Let po E L 2 (T) be the value of p at xq. 
Let us assume that B and Q are decoupled from the space variable i.e. they are such that: 



S(x - x Q )Bg('y,p)dx = Bg(~f,p ). 



X 



It is the case for the clarinet where the air flow Q at xq depends only on the pressure in the player's 
mouth (7), the internal pressure p at xq and on the size of the opening at this point. Let us further 
assume that A and Q possess the properties allowing to write: 



p(x,t) = / p(x,t-T;Z)6(Z-x ) BG(~f,p)\ i T d£di 



XxT 



J^p(x,t-T;x )Bg( 7 ,p )\ T dT. (2) 



IT 

As a consequence of the above equation, there exists an operator Z : L 2 (T) — > X such that: 

p = Zg(~f 1 p ). (3) I eqPZG 

Notice that the action of Z implies ac^gryolution with respect to time, which becomes a multiplication 
by Fourier transform. Expressing §3| at x — xq, one can introduce the operator Z = Z\ xo , where 
Z : L 2 (T) — > L 2 (T), and the study is hence reduced to the investigation of the oscillating solutions 
belonging to L 2 (T) of: 

p = Zg(^,p ). (4) eqreel 

Equation |T]) is a formal expression whereas (JSJ) is rigorously defined in the following and is the starting 
point of the study. The woodwind instruments have been modelled by physicists in the framework of 
control theory where the equations are written in term of transfer functions. These latter quantities are 
the Laplace transform of some differential equations, and differential operators bgcom^ multiplication 
operators under thi^ transformation. The Fourier transform of the operator Z in @ is a multiplicative 
operator and @ follows the control theory formalism. Thus, the present study is oriented toward con- 
trol theory, however, with appropiate assumptions, it also applies to the usual description with partial 
differential equations. 

The nonlinear vector valued function Q : (R x L 2 (T)) — > L 2 (T) has the particular shape: 

00 

G(l,Po) = Lp + ^ R n(l,Po), (5) I eqdep2 

n=2 

where the linear operator L depends on 7 and each R n corresponds to a monomial of degree n in po. 
More precisely, for each n > 2, let Z n< \, Z n> 2, ■ • • Z n , n be bounded linear operators from L 2 (T) to itself. 
Each operator may depend on 7. It is assumed that the nonlinear terms can be written as 

Rn{l,Po) = (Z 7h lPo)(Z rh2 po) ■ ■ ■ (Z ntV ,Po). (6) 

This expression of the nonlmearity allows to cover the case of woodwind musical instruments as it will 

r t Tw oodwind 

be shown in section |[7] as well as different types of delay differential equations. 
Remark 2.1. Let some dynamical system be described by the standard notation: 

^jjjT = L 2 p + 9(7, Po) (7) I eqstd 

IGH 

where L2 is theJinear part and g the nonlinear mapping. See \ \'ZU\ for examples of assumptions on 
theses objects, 10 for the infinite dimensional case, or 1 (1 1] for the case of the Hopf bifurcation. The 
connection with this study is made by writing formally Z _1 (l — ZL) = L2 — d/dt. By Fourier transform 
the time derivative operator becomes the product by some complex number iui, lu E M. Looking for a 
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pure imaginary eigenvalue of L2 is the same problem as searching for a zero eigenvalue of Z~ x (\ — ZL). 
In order to prove the existence of a Hopf bifurcation, the norm of the resolvent (L^ ± must be 

bounded by some constant/ |w| when the spectral parameter ±iuj tends to complex infinity (see l|3]), a 
property which does not arise in hyperbolic PDEs. Here the inverse of Z~ x (\ — ZL) is only required to 
be bounded for each 7 in a neighborhood of the bifurcation point, except at this critical point. However, 
this decrease of the norm of the resolvent is replaced here by the assumption of a twice continuously 
diffcrentiable periodic solution. 

Let L°° (T) be the space of bounded functions on T and let B e C L°° (T) be the ball of radius e > 
centered at zero. Let C 2 (T) be the set of functions twice continuously differentiable with respect to t € T 
and define 

A = 6 £ nC 2 (T). 

The solution of (|[]) under investigation satisfies the following hypothesis: 

HI For all e G (0, 1) there exists a non-empty, open, simply connected set Qi Cl, where for all 7 G fii 
the system has a unique non trivial real-valued periodic solution belonging to A, continuous with 
respect to 7. Its period will be denoted by t = t(j). There exists a 70 G R, with 70 G fl the closure 
of fix, such that the periodic solution tends to zero when 7 tends to 70. The point 70 will be called 
the critical point or the bifurcation point. 

As said in the introduction, and stated in HI, the regularity of the solution is an hypothesis. Nev- 
ertheless, several publications have shown regularity properties of the solutions of nonlinear hyperbolic 
PDEs. For instance IT2J shows that a twice continuously differentiable solution with respect to time exists 
for a model of the wave equation with damping and source term, this solution being a distribution with 
respect to the space variable. In \ [16\ , for a slightly different model, the solution is twice continuously 
differentiable with respect to time and, for all time fixed in a bounded interval, it belongs to a Sobolev 
space. The periodic solution^ of the nonlinear wave equation are even more regular for some class of 
nonlinear terms as shown in IT4J and references therein. To the knowledge of the author, no such results 
have been sett ^edjba the case of a nonlinear function localized at one point in space. However, the model 
described by £§q Skfe^ ^ represents J^jp ^ave equation in a tube with damping, has been investigated 
theoretically |T67 1 Y J and numerically |T87"E9| by physicists, giving results compatible with this regularity 
hypothesis, and motivating the present approach. 

Let us introduce the notation 

J 7 = Z~ x (l - ZL), 

and assume: 

H2 I) J 7 is a closed operator acting in L 2 (T). It is diagonalizable and its eigenvectors consist of the 
vectors of the Fourier basis of L 2 {T). 

2) Let Q e , Q_ e be the orthogonal projections associated to the vectors ip e , ip- e , respectively, of 
the Fourier basis. J 7 is such that 

QeJ-yQe = (3 e (j)Qe, Q- e J~{Q-e = (3-e(l)Q-e 

where (3e(l), (3-e(l) are two complex conjugate eigenvalues of J 7 such that e {lo) — 0. More- 
over, for all 7 € fi, the operator QJ 7 Q, where Q = (1 — (Q e + Q- e )), possesses a bounded 
inverse. 

On the nonlinear part, the hypotheses are: 

H3 For all 7 e il and n,m £ N, n > 2, m € [l,n], Z n , m (-f) : L 2 (T) — ► L 2 (T) is a bounded linear 
operator. Moreover, as for J 7 , the Fourier basis is its eigenbasis. For all n > 2, 7 G £1 and p e A, 
by the Holder inequality, the quantity i? n (7,po) belongs to L 2 (T), since the {Z niTn } are bounded 
and ||po|| 00 < e, where || • ||oo denotes the norm associated to L°°(T). Eventually, it is supposed 
that for all 7 G D, and po € A, the series of nonlinear terms is convergent : 

M 

lim V R n (j, P() ) = R(7,p ) G L 2 (T). 

M — s-oo * — * 
n=2 
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Remark 2.2. If for all 7 £ f2, Q is an analytic function in its second variable, then it satisfies the above 
hypothesis: all the Z n ^ m are identity operators except when m = 1, where Z n ^\ is the term of the Taylor 
expansion Z n _\ = Q^> (7, 0)/n!, with C/^™" 1 denoting the nth derivative. 

Since the solution is periodic in time, it is natural to work in the Fourier space. Let S T be the circle of 
perimeter r, then T = S T . Let us define £ 2 , the Hilbert space of sequences which are square summable. 
Let T denotes the Fourier transform operator from L 2 (T) to I 2 . The Fourier transform of po will be 
denoted p~o = !Fpo £ i 2 or, using the physicists's convention, by a capital letter P = §q. For all q £ Z 
and Lo q — -^r-, the Fourier transform is defined by : 

P q =p(cj q ) = (T Pa )(Lu q ) = - [ pofle-Wdt, 

and the inverse transform at time t £ S T is: 

{T- 1 P){t) = ^P q e lu) « t . 

Using integration by part in the expression of the Fourier transform, the hypothesis of twice continuously 
diffcrcntiablc solutions of HI implies the existence of ao > such that: 



Vg ? 0, \P q \ < 



M 2 



(8) 



Pdecr 



Moreover, from HI and the definition of the Fourier transform, the small amplitude condition gives for 
all 7 £ r2 ; 

Sup|F 9 |<£. (9) I boundPql 



3 Expression in the frequency domain 

Many physical systems having an oscillating behavior, for instance in mechanics, acoustics or electrical 
engineering are described by relations taking place in the frequency domain. This is the formalism of 
transfer functions. Moreover, investigating oscillations turns out to be extremely efficient when using a 
frequency domain approach. In finite dimension, the harmonic balance method or the describing function 
method are such tecl nikm c^. In infinite dimension, these techniques have not been proved to hold. This 
section, with lemma lb. 11 shows how the infinite dimensional system under study can be written in the 
frequency domain. It is the first step toward its analysis using frequency domain techniques. 

The key point is to express the nonlinear relationship in the Fourier basiSpJt c will become a combination 
of convolutions. Let g,h £ £ 2 be such that their harmonics decrease like in (0), m consequence they belong 
to e 1 , the space of summable sequences. Thus, the convolution g * h defined at q £ Z by : 

{g*h)(q) = ^2 K-n9n, 



is an I 1 function by Young's inequality, see e.g. I"2T, Sec. IX. 4]. For each Z n<m defined in H3, let us 
introduce the notation: 

Zn.m — 3~ Z n m J~ 

From H3, since each Z n:Tn is diagonal in the Fourier domain, one has for the value of P associated to the 
vector i\) q of the Fourier basis: 

{J 7 Z n>m p)(q) = Z n:Tn (q)P q 

where Z n ^ m [q) is a complex number. 

More generally, let p, u±,U2, • • • , ujv-i be N > 2 real periodic functions depending on time t £ S T . 
Let these functions describe the state of some physical system under study. Let P, , , • • • , t/^ -1 ) 
be the Fourier transforms of p, u±,U2, • • • , wjy— 1 respectively. Suppose they are related together through 
N — 1 linear relations, in the frequency domain, involving transfer functions: 



t/(") = Z m P, 



(10) 



UZP1 
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where the {Z m } are bounded operators. All may depend on the bifurcation parameter. Let us suppose 
there is also one nonlinear relation, often expressed in the time domain: 



Af(j,p,U!, it 2 , • • • ,ujv-i) = 0. 



(11) 



If this relationship involves polynomial terms between the variables, such as p 2 , pu\ or u\u 2 , the following 
lemma can be used to express this relation in the frequency domain, with convolution products. For 
example, a cubic term involving three functions p, u\ — Z\p, u 2 — Z 2 p, with p satisfying HI and Z\, Z 2 
satisfying H3 would give : 

pwrui = p*u{ *u 2 ~ = p* {Zip) * (Z 2 p). 

And for all q £ Z: 

jp* (Zip) * (Z 2 p)| (U)g) = ^2 X! Z ^ UJ k 1 - UJ k2 )Z 2 {uk 2 )P q -k 1 Pk 1 -k 2 Pk 2 - 

fciGZ fe 2 ez 

I nonlinex 

The whole rdation can be stated under a set of harmonic balance equations, in the frequency domain, 
given by flop . 

Lemma 3.1. Suppose pq £ A and P denotes its Fourier transform. Under the hypotheses H2, H3. for 
all 7 € fi, all n £ N, n > 2, 

TR n {-i,p Q ) = S(g n ,n,P), 
where S(g n ,n, P) £ £ . For all q £ Z the elements of S(g n ,n, P) are defined by : 



Sq(9n,n,P) = X] Z\2 "* X! 9n(q,h,--- ,k n -l)Pg-k 1 Pk x -k 2 ---Pk n - 1 , 

where for all n, g n is a complex-valued function of n variables depending on combinations of Z n ^ r 

g n (q, fci, • • • , kn-x) = Z n> x(q - ki)Z n)2 (ki - k 2 ) ■ ■ ■ Z„, n (fe„_i). 
Moreover, the infinite sum gives: 

oo oo 

FYsM'YiPo) = S (9n,n, P)££ 2 . 

11=2 n=2 
Jeoreel 

In consequence, solving the system described by (Q) is equivalent to solving: 

oo 

T 1 P = Y J S{g n ,n,P). 



(12) 



nonlinex 



defS 



(13) eqfreq 



Notice that although not specified, the convolution S depends on 7. 
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Proof of lemma W. it -Let G^ 2 \ • • • , G'" +1 ^ £ i 1 , n £ N, n > I, then the n-product of convolution 
G« * GP) * • ■ ■ * G(" +1 ) is an t 1 function, by Young inequality. Suppose g,h £ L 2 (T) such that their 
respective Fourier transform G, H are i 1 functions. Then Fubini's theorem can be used to get: 



q£Zk£Z q&Zk&Z fc£Z g£Z 

= Y,Ht)Gke iUkt = h{t)g{t). (14) 

fcGZ 

Notice that the functions h, g are continuous by the Riemann-Lebesgue lemma. Reciproralh/, for h, 
g £ L 2 (T) twice continuously differenti^lej i? and G are i 1 functions (see explanation of (0) j and this 
equality holds. In consequence, Eq. (|IL-4|) snows the correspondence between multiplication in the time 
domain and convolution in the frequency domain of a finite number of products, in the case of twice 
continuously differentiable functions on the compact domain T. Since H3 holds, all the {Z n ^ m } are 
bounded and for all n > 2, g n is bounded. The bound may depend on n. Using Holder inequality, one 
gets: 



\\Rn(i,p)\\mT) < \\z n ,i 



z. n 



(Ml 



IHU 2 (t), 



mulconv 
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which shows that R n ( / y,p) belongs to L 2 (T). The symbol || • || denotes the operator norm. Suppose 
S(l,n,P) £ £\ then 

5(1, n + 1, P) = P * 5(1, n, P) e f 
and since 5(1, 2, P) € i 1 , by induction, for all n > 2, 5(1, n, P) G Then, one has: 

\\S(g n ,n,P)\\ e i < sup\g n \ ■ ||5(1, n, P)\\ t t. 



In consequence, (|ll4|l snows that for n > 2 hxed, S(g n ,n, P) is the Fourier transform of R n (j,p). Let us 
study the limit n — > oo. Since the Fourier transform is bounded, and even unitary one has 



S(P)=FR(>y,p)e 



and the quantity: 



M 



S(P)-J2S(9n,n,P) 



< 



t 2 



M 



L 2 (T) 



tends to zero as M tends to infinity from H3. 



□ 



threshold 
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4 Preliminary results and properties at threshold 

This section demonstrates how the method based on the frequency domain approach can be used to 
investigate a bifurcation in a nonlinear system. It shows how to handle the hypothesis and the Fourier 
harmonics in order to control infinite sums and convolutions. The first result stated in the following 
theorem concerns the determination of the birth point of the oscillation and the need of a zero eigenvalue 
for the linear part. Notice that it could have been directly given by means of the Hartman-Grobman 
theorem, however, it was not the case so as to demonstrate how the frequency approach works. The 
second result of the theorem gives an estimation ^ t he harmon^s^and shows the importance of the first 
harmonic. This serves as an introduction to Th. lb. 1 1 of section Q which will go further in this estimation. 
Finally, this section settles many important intermediate results which will be used in the other parts of 
the study. 

Theorem 4.1. Suppose HI, H2 and H3 hold. Then the bifurcation point 70 where the oscillating 
solution starts and the frequency u e of oscillation at this point are given by the characteristic equation: 

A, (70) = 0. 

Moreover, for all j £ Q where e is small enough, \P e \ — e and there exists cq > such that for all 
qeZ\{e,~e} , \P q \ <c e 2 . 

Remark 4.2. The q^Sti^ff e i £ an d 7 are dependent of each other and this dependence will be made 
explicit in the section |dJ 

Remark 4.3. In the case of the Hopf bifurcation, the condition of two purely imaginary eigenvalues of 
the Jacobian operator at the bifurcation determines 70 and the oscillating frequency lu c at this point. In 
the present case, the result is equivalent with the condition on two complex conjugate eigenvalues of J 7 
to be zero. 

Before stating the proof of this theorem, several lemmas giving bounds on the different convolution 
products encountered are presented. Let us start by showing a property of the convolution which will be 
crucial in the following. 

Lemma 4.4. Let ai, 02 be positive constants, X € (0, 1) and 9%, 62, oti, c<2 be strictly positive parameters. 
Let f,g:'L'^ ( C be two bounded decreasing functions such that 



VqeZ, \f(q)\<X e \ \ 9 (q)\<\ 6 \ 



and Vg G Z* 



1/(9)1 < 



Q2 

,|1+C*2 ' 



Then. 
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i) there exists a positive continuous function c\ : (0, 1] — > M + such that 

Vg G Z, |(/ *.g)(g)| < mm{ a2Cl (a 2 ) • A ri ,a 1 c 1 (a 1 ) • A r2 }, 

with 

n = 1 + — — , r 2 = t> 2 



2 1 + a 2 ' 2 1 + ai' 

iij £/iere exists a small a > and c 2 = c 2 (a, ot\, a 2 ) > swc/i iftai /or A smafi enough, Vg G Z* 

|(/*<?)(?)| ^-^A^A 8 ^-^, 

with r = 1+a ^ +a2 - a(2 + ai + a 2 ). 

As a consequence, for any number M > 2 of functions /i, / 2 • • • , /m satisfying the hypothesis for 
some and {ai} il i G [1, M], the multiple convolution /i * / 2 * • • • * / M is a bounded function with 

bound depending on a positive power of A. 

Proof. Let us introduce the notation f q = /(g) for / : Z — > C. Let /, g be two functions satisfying the 
hypothesis. Let oi, <r 2 G (0, 1] such that (1 + ai)(l — o"i) + (1 + a 2 )(l — er 2 ) > 1. For all q G Z, 



< r > ue 1 +e 2 ) , \0i<xi+02<T2 i-^i i-<72 1 

" A +A 0l ° 2 2, |fl_W(l+ai)(l-«ri)|jfe|(l+a a )(l-«r 2 ) 



< a{- ai al- a2 c(o- u 0-2, a 1 ,a 2 )\ 9iai+e ^\ 

where c(a\, <r 2 , ai, a 2 ) > tends to infinity when o"i, cr 2 , ai, a 2 are such that [(1 + ai)(l — o-\) + (1 + 
£* 2 )(1 — cr 2 )] tends to unity. The last step in the above calculus is obtained using Holder's inequality. Fix 
a i = 1, then the convergence criterion implies a 2 < 1 . If cr 2 is fixed to the particular value 

1 Q! 2 . . I ; [ 

<T 2 = -— , (15) fixsigma 

21 + a 2 I 1 

then the constant c depends only on a 2 and there exists a constant Ci(a 2 ) > such that the convolution 
is bounded by: 

\y2f q -k9k\ < a 2 ci(a 2 )X ri , n = 6\ + 6 2 \ ^ . 

feGZ 2 1 + a2 

And similarly for cr 2 = 1, one obtains the constant aici(ai) and the value r 2 = # 2 + ^i ^ ■ 
For ii), for g G Z* , a G (0, 1 "^"^ ) the following estimate can be written: 

where 

5i = (l + ai)(l-a), i = 1,2. 
Fix (3 — min{<5i, <5 2 }, for any er G (0, <5l+ 2 < ^~ 1 ) one has the bound: 

v 1 1 V 1 1 



< c(a, a, ai , a 2 ) ^sup^ ^ _ 1__ , with c(a, a, a lt a 2 ) > 0. 
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The calculation of the minimum over k of h q (k) = \q — k\\k\ will finish the setting of a bound for the 
above quantity. For k ^ 0, k ^ q, the derivative h! q is h! q (k) = — sgn(q — k)\k\ + sgn(fc)|q — k\. It has the 
following sign for q > 0: 



k 


-oo 


q/2 q 


oo 


K 


- 1 H 


- - | - 


- 



If q < 0, taking (— k) instead of k leads to the same result. The function h q has a minimum at k — and 
k = q (where it reaches zero) and a maximum at k = q/2. Hence, for k,q € Z, k ^= 0, k =/= q and \q\ > 1, 



the minimum of the function is: 



min ft g (fe) = \\q\ - 1|. 

kGZ,* ,k^q 

Remark that for q — or ±1, the minimum is 1 as k ^ q and k ^ 0. Moreover, for all g such that \q\ > 1, 
one has \q\/2 < \\q\ — 1| and 



E 



i 



k£Z*,k^q 



\q- k\ 5 i\k\ s * 



< 



2 l3a c(a, a, ai, 02) 



1/9(7 



<5l+A" 2 -l _ l+ai+q 2 _ 2+ai+q 
2 2 u 2 

exists a small 77 > and a positive value 02(77, a, ai, 0:2) such that taking (3a = S — 77 gives 



For the case q = ±1, the bound is simply c(a, a, a 1 ,a 2 ). Let S = °i±2a=A = ++"i+"2 - a 2+a ; + " 2 , there 



1 e < 1^ (( aifl 2) i -A^^)" C ~ 2( ,,a,a 1 , a2 ) + |g|(1+ ;; A !; )/2+ , + k|(1+ :; A ! 2 1)/2+ , 

This holds for the case g = ±1 as well. Since a is small, (a^) 1 "" < 2aia2 and A 1 < A ia . The parame- 
ters a and 77 are arbitrary small so one can choose a{2 + u\ + ct2)/2 = 77. Let di2 = max{aiA e ' l0! , 02 A 6 * 2 ™}. 
For A small enough, the first term of the above estimate is the largest one and there exists ci(a, ot\, a%) > 
such that 

E, I < di2C 2 (a, ai, ag) 

Jq-k9k\ _ 1 |(i +Qi+Q2 )/2-q(2+ Qi +Q2) ' 



□ 

The next lemma states upper bounds on the value of the function S defined in 
Lemma 4.5. Suppose there exists ao > 0, and e > such that for all q G Z, 

»; 1^1 < e. 

4l J /or allqeZ*, \P q \ < 
Then there exists a small a > suc/i i/ia£ /or a(In6N,i!>2, <;£Z and e small enough: 

\S q (l,n,P)\ < e i+("-i)(i-°). (16) 

JlemmeSl 

Proof. Let us first investigate the case n = 2. The hypotheses i) and ii) allow to apply lemma B. 41 with 
0\ = 62 = 1, ol\ = 1, ot2 = 1, ctx = ci2 = cto, A = e. This gives: 

VqeZ, |S,(l,2,P)| = |(P*P)(g)| <a Ci-e 5 / 4 , 

1 



VoeZ*, |S,(l,2,P)|<C 2 (a)age 



where Ci > 0, a > is small and C 2 (a) > 0. For e small enough, the quantity C2{a)a^e a is lower than 
one and aoC\e a < 1. So that the previous relations become: 



VgeZ, \S q (l,2,P)\ = \(P*P)(q)\ < e 5/4 " 



VoGZ*, |5,(1,2,P)|< 



1 



|i+l-4a' 



(17) 
(18) 



estimS 








estimSSS 



9 



Let us remind that the functions 5(1, i, P), i > 3, can be written for all q G Z 



lemmeosc 



Using lemma 11.41 with this formula and (|llV|) . l|ll8]) . this yields the following bounds for 5(1,3, P) 



5(l,i,P) = P*5(l,i- 1,P). (19) 
nd (|lt V|) . l|ll8]) . this yields the following bounc 

Vq G Z, 3, P)\ < aoCi • £ 5/4-«+l/4 < £ l+2(l/4-a) ) 

VgeZ*, |5 9 (1,3,P)| < C 2 (a)a e a • ^- < ' 



Srecur 



t q\ u 3 - \q\ u 3 

where u$ — 1 + 1/4 — a(5 + 1/2) + 4a 2 . Similarly, by induction, each quantity 5(1, n, P) is bounded by: 

V?eZ, |5 g (l,n,P)|< £ *», (20) 
VaGZ*, |5 g (l,n,P)|<^, (21) 

with v n = 1 + (n — 1) (1/4 — a) and u n = 1/2 + u n _i/2 — a(2 + u n _i). Remark that for all n € N u n > 
and 

2a 

hm m„ = 1 . 

n— too 1 — 2a 

□ 

Jsaf req 

One can now give a first estimate on the harmonics P q of the solution of 

Lemma 4.6. Suppose HI. H3 /io/d. For e small enough, relation (|I13|) implies the following properties: 

i) if hypothesis 1) of H2 holds and the operator J 7 is invertible for all 7 G fii, J 7 mus£ noi 6e 
invertible at 70 . 

izj z/ H2 holds, then for e G Z* defined in this hypothesis, P e ^ and e can be chosen to be e = \P e \- 
Moreover, there exists a constant cq > such that, for all q G Z. q ^ e, \P q \ < Co£ 2 . 

zizj i/ H2 holds, for all (3 G Z smc/i i/iai /3 ^ { me } m eZ' -^/9 = 0- 

Proof. Suppose J 7 has a bounded inverse for all 7 G fi, then the complex- valued function J 7 is boundec 
Since H3 holds, the functic 
gives for all q G Z, M > 2: 



estimSn 



estimSnq 



freq 

Since H3 holds, the functions {g n } are bounded. Let us write g n for the supremum of |<j>„|. Equation 



1 M 
|P g |<sup|J 7 (k)\ X sup <?„ x I V 5 9 (l,n,P)| + 

ft ne[2,M] n=2 

oc 

+ sup|J 7 _1 (fc)| x I S q (g n ,n,P)\. (22) 

fc n=M 

J lammetf 

From H3 and lcmma lB.il tlie second term of the previous expression tends to zero as M tends to infinity. 
Let us choose M large enough such that 

00 M 

I V S g (g n ,n,P)\ < sup g n x | V 5 9 (l,n,P)|, 

n=M ne[2,M] n=2 

then there exists a constant ^4 > such that 

M 

\P v \<Ax\J2S q (l,n,P)\. (23) 

n=2 

JlammeSnl 

Since HI holds, the hypotheses of lcmma B.til are satisfied and for e small enough: 

Vn> 2,qe Z,|5 9 (l,?i,P)| < £ l+(»-l)(l/4-a) ) (24) 



estimJP 



lboundS 
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where a is small. Thus for all g€Z: 

M 



M-2 



\P q \ <AeY,£ (n ~ m/4 ~ a) = Ae 1+1 / 4 - a £ £ ™d/4^) ; 



n=2 



m=0 



and for e small enough such that Ae a X^m=o £ m ^^ 4 Q ' < 1: 



l^l< £ 1+e , 



= (1/4 -2a). 



(25) 



(26) 



boundP 



eSnl 



For i), using this new bound as hypothesis for lemma B~bi ( replacing e by s 1+9 ), gives for all q £ Z a new, 
stronger, bound on P q , i.e.: 



|P,| <e {1+6) 



(27) 



boundP2 



Since > 0, doing this reasoning again, an infinite number of times, will lead to P q — 0, for all q. 

For ii), let s be the set s = {e, — e, q — e,q + e}, where e is the integer defined in H2 and geZ. t ^gi'e iln1p 
J 7 is not y^vertible for all 7 G f2. One can replace J" 1 by (QPyQ) -1 which is bounded and estimates tflZ'd)) 
and (p6|l liold here, for all q ^ ±e: 



|P,| <e 1+e . 



The function 5(1, 2, P) at q S Z can be written as: 

2, P) =2P 9 _ e P e + 2P g+e P_ e + J] P g _iP. 



(28) 



(29) 



boundP22 



devS 



dP22 



Using paf) . trie first two terms are bounded by a quantity proportional to e 2+s in the case where q 7^ ±2e 
and <? 7^ and e 2 otherwise. The sum over n contains only temis^ljounded by e 1+e so has a bound 
proportional to e to the power (1 + 9) 2 , consequence of lemma ft. SI Similarly, for q € Z and n > 3: 



5,(1, n, P) = (P * 5(1, n - 1, P))(<?) =P e S q - e (l, n - 1, P) + P_ e 5 9+e (l,n - 1, P)+ 

+ P 9 _ e 5 e (l, n-l,P) + P g+e 5_ e (l, n - 1, P)- 

+ 5^P,_ i 5 i (l,n-l J P). 



(30) 



its 



JdeyS JlampeSnl 

This is just %L\)\ where the Pj have been replaced by 5j(l, n — 1,P). So, similarly, using lemma f4.a|) . tor 
e small enough, it implies: 



Vn>2, geZ, |5 g (l,n,P)| < £ (1+e 



(31) 



imJP 



The qua ntity jgj^jp? is invcrtiblc, from H2, then for all q 7^ ±e one has the estimate (|l^3|l and for e small 
enough, if) gives: 



IPI <£ (1+9)2 



(32) 



Similarly to the proof of i), this procedure can be repeated s^v e^al times and gives each time a larger 
power of e, until the first two terms of (|Li!j)) and four terms of (|bu|) become the largest terms. This implies 
that for all q € Z, q 7^ ±e, \P q \ < |P e |. If P e = 0, then no oscillation exists. If P e 7^ 0, the largest term 
P e can be chosen to be \P e \ = e and there exists cq > such that for all q € Z, q 7^ ±e 



IPJ <c £ 2 



(33) 



The above estimate is given and limited by the case q = ±2e. 

For Hi), let us define the set s e = {me} me2 ,. If /3 ^ s e , then for each n e Z fixed, at least one of the 
quantities (f3 — n) or n does not belong to s e . Suppose for all n £ s e . 



devS2 



bounds all 



boundPq 



IPJ <e\ a> 1. 



(34) 



steplP 
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The first two terms of (pup are bounded by e a+ and from i) of lemma E. 41 one can write: 

with a > small. It is similar for (pup . Following the procedure used in the proof of ii), one will find 
that a is increasing at each step. Eventually, a tends to infinity and implies for e small enough: 

P =O. 

□ 

Proof of Theorem \\4.1\ The proof is a direct application of the results given in lemma lB.bl INotice first that 
in order to have a solution satisfying HI, the point i) of this lemma requires that J 7 possesses at least 
one zero eigenvalue at 70 G f2, showing the necessity of a hypothesis such as H2. As a consequence, the 
bifurcation point is the point for which an eigenvalue of J 7 reaches zero. The point ii) of the same lemma 
states that the Fourier series of the appearing oscillation must include a non zero value for P e , where e is 
given in H2. The point Hi) implies that P e is the term of the Fourier series associated to the frequency 
of oscillation, since only harmonics associated with frequencies multiplc,,of e can be different from zero. 
Eventually, the bound on the Fourier components given in ii) of lemma B.bl finishes the proof. □ 



decrease 



Thl 



5 Decrease of the harmonics 

This section concerns lhe estimation of the amplitude of the Fourier components of the solution. It is 
dedicated to theorem lb. II and its proof. 

Before stating the result, notice that the index e of H2 can be chosen to be e = 1 without loss of 
generality. This is justified by ii) and Hi) of lemma 16. bl wnich asserts that the fundamental frequency is 
given by u> e and there exists only harmonics associated to multiple of this frequency. 

Theorem 5.1. Under the hypothesis of Theorem \\4-l[ for all 7 G ft where .j-fjj is small enough, there 



exists k = fc(|Pi|) > 2 such that the oscillating solution of the system 
properties: for all ggZ*, 



\P q \<c q \Pi, 
1^,1 < c*|Pl| 



\q\ 



for all \q\ < k 
for all \q\ > k, 



JraJ vossesses the following 



(35) 



estimPk 



where c q ,Ck are positive constants. The rank k increases when \Pi\ 
I Pi I tends to zero. 



decreases and tends to infinity when 



Remark 5.2. If some of the transfer functions {Zij} change sign or decrease fast enough when their 
variable tends to infinity, the decrease of the harmonics amplitude with frequency may be faster than the 
estimation given in the above theorem. 



Proof of theorem The reasoning of the proof is similar to the one used in the proof of lemma Ikt.bl 
The operator Q = (1 — Q\ — Q-i) commute with the Fourier transform since Qi, Q_i are associated to 
vectors of the Fourier basis. Then the operator QJ 7 Q has a bounded inverse for all 7 G O. Since H3 
holds, the functions {g n } are bounded. Let us write g n for the supremum of \g n \. Equation 3p gives for 
all q G Z\{1,-1}, M > 2: 



A I 



\P q \ <sup|(<3J 7 Q) ^fc)! x sup g n x 

k n£[2,M] 



J2S q (l,n,P)\ + 



- sup I (QJ^Q)- 1 ^) I x I J2 S q (g n ,n,P)\. 



(36) 



netf 



J lgmme 

From H3 and lcmma lB.il the second term of the previous expression tends to zero as M tends to infinity. 
Let us choose M large enough such that 



M 



\ S q (g n , n, P)\ < sup 

'J n 

ne[2,M] „ , 



n=M 
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then there exists a constant A > such that 



M 



\P q \<Ax\J2S q (l,n,P)\ Vflf^il. 



(37) 



n=2 



Let s be the set s = {0, 1,-1,2, -2,q- l,q+ l,q- 2,q + 2}. The function 5(1, 2, P) atgeZ can be 
written as: 



5,(1, 2, P) =2P,_ 1 P 1 + 2P g+1 P_! + 2P 9 _ 2 P 2 + 2P, +2 P_ 2 + ]T P 9 -„P„ 



(38) 



™ez 



From Th. lE II tor all |g| > 3 the first four terms of (|i^J) are bounded by a quantity proportional to e 3 . 
Likewise, but for all q G Z, each texmoOhe sum over n contains two quantities, each one is bounded by 
cq£ 2 for some cq > 0, then lemma 16.41 gives the estimate: 



^|P 9 _„P„|<(^) 2+1 



/2-c 



(39) 



n£Z 



with a > small. As a consequence, for |g| > 3, and for e small enough, |5,(l,2,Pl| S^£ 2+e where 
9 = 1/2 — 2a. In the case where q — 0, q = ±1 and q = ±2, some of the first terms of (|3&|) are bounded 
by e 2 , so for q G {0,±1,±2}, there exists c 2 , 2 > such that |5 g (l,2, P)| < c 2j2 e 2 . Let us now give a 
bound on the 5(1, n, P) by induction. Assume that at rank n — 1 > 2 the following statements are true: 

• forge {0,±1,±2}, \S q (l, n - 1, P)| < Cie 2 , 

. forg, |g| > 3, |5 ? (1, n - 1, P)| < e 2+e , 

for some Ci > and 6 — 1/2 — 2a. Then at rank n, for g G Z : 

5,(1, n, P) = (P * 5(1, n - 1, P))(g) ^U^ 1 + ^ P,-^!, n - 1, P). 



(40) 



Its 



where, for s 2 = {0, ±1,±2}, 



K' 1 = P kS q -k(l, n - 1, P) + P,_ fc 5 fc (l, n - 1, P). 

fc£s 2 

Firstly, this latter e^pr,gssion is bounded by a quarrtitv^roportional to e 3 for all g. Secondly, for all g G Z 
the sum over i of (Ifauj) is bounded, from lemma 16.41 bv: 



E 



P 



8-2 



CO 



5j(l,n- 1,P) 



< 2+8 + 1/2-a 



Then for e small enough, at rank n: 



VgGZ, |5,(l,n,P)| <e 



2+26 < £ 2+9 



imJP2 



and by induction this relation holds for all n > 2. In consequence (|bV|) implies that there exists a constant 
D > 0, such that for all |g| > 3: 



estimJP2 



devSp 



devSS 



\P q \ <De 



2+6 



(41) boundPq3 



Similarly to the proof of ii) of lemma lb 61 this proceclure can be repeated several times and gives each 
time a larger power of e, until the first terms of (|l3§|) and some of the terms of U q become the largest 
terms. The bound on these latter terms is here proportional to e 3 . In conclusion, there exist P/ 2 ,P/3 > 
such that one has: 



forgGs 2 , \P q \ <D 2 e 2 , 
forg£s 2 , |P,| <P 3 e 3 . 

For the general case, suppose there exists a rank r > 3 and a function / : Z 



(42) 



such that 
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. for all q, \q\<r, \P q \ < f(\q\)e^, 

• for all q, \q\>r, \P q \ < f(r)e r . 

Let us show it is true at rank r + 1 under appropriate conditions. The steps are the same as for the case 
r = 3 treated above. In the following, only the case q positive is treated as the case q negative will lead 
to the same estimates. For g g [2, r + 1], let us introduce the set s q = Z fl [1, q — 1] and the function / 
defined by : 

Vq>2, f(q) = J2f( ( l-i)m, 7(1) = 1. (43) 

Notice that the function / is monotone increasing on N. Let us introduce for q > 2 and n > 3 

u l = E P i-^ U 2 = E P 9-iSi(l, n - 1, P). 

i&s q i£s q 

Instead of (|b»|) . one has to estimate the quantity: 

S q (l,2,P) = U% + J2 P <i-i P i- ( 44 ) 

For q < r, the quantity U 2 is bounded by: 



\U 2 \ < £ ^/(g- »)/(<) = /(«)** 



eSl 



5 g (l,2,P)=^ a + 5^P,_ i P i . (46) 

The quantity V q 2 contains a small number of terms bounded by e r+l , this number is independent of q, 
and the others are bounded by e r+2 or lower bounds. So for e small enough, there exists a g?o > such 
that: 

\V q 2 \ < d f(r)e r +K 

For the ^um ^oyer i of (|m>)) . each term contains two quantities where at least one is bounded by f(r)£ r . 
Lemma [Oj gives: 

lE P ?-^l < ci/W 2 e r+1 / 4 -", 

where 77 > is small. Let us introduce 

/ m (r + l) = max{/ 2 (r),7(r + l)}, 

then there exits c?2,2 > such that for all q > r + 1: 

|5 g (l,2,P)| < d2,2/ m (r + l)e r+1 / 4 -". (47) 

Next, one has to treat the case of multiple convolutions. Similarly to the case r = 3, suppose at rank 
n — 1 > 2 there exist a function h n -x : Z — > M+ such that: 



Usumi 



Jllaiuiii 11/1 . ijjimiu^^ ^ 

and the sum over i of OFgrTs bounded by e 9+Vi-« j for some small a > 0, from lemma (4.4. Hence, for t 

small enough, there exists c?2,i > such that: 

Vq<r, \S q {l,2,P)\<d 2 J{q)e". (45) [^17 

For q > r + 1, let us introduce the set t r — Z n [— r, r] and for n > 3 

y 9 2 = E p 9-* p *> y <r = E p i-i s i(^ n - i. p )- 

Instead of (p&|) . one has to estimate the quantity: 



SVsum 



qsupr 
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. forge [-r,r], \S q (l, n - 1, P)| < d„_i,i/i„_i(|g|)el9l, 
• forg, |g| >r+l, |5,(1, n - 1, P)| < d n _ 1)2 /i„_x(r + lK+ e , 
for some d n -i,2 > and € (0, 1/4). Then, at rank n, for |g| < r : 

S g (l, n, P) = (P * S(l, n-1, P))(g) =E7?- 1 + ]T P,_<Si(l, n - 1, P), (48) 



devSSg 



i<£s„ 



and for |g| > r + 1 : 

S q (l,n,P) = (P*S(l,n-l,P))(q)=V^- 1 + J2 p i-M^ n -^ p )- ( 49 ) | devSS g 2 



is£t r 



It is similar to the previous case, where n = 2, provided £ is small enough. More precisely, it is required 
that for all n & [3, M + 1], one has e such that 

/i„_i(r + l)e < 1. 

Then the quantity U™ is bounded by 

\U™\ < K{q)e\ where h n (q) = f(v ~ 



The term V^Ms bounded by a quantity proportional to h(r)e r+1 for e small enough, the sum over i of (IBS' 



(resp. can be shown to be bounded by a quantity proportional to el <? l+ e (re sj^ , g^^)- This finishes 

the estimates on the multiple convolutions. The last step of the proof is to use (pVj) and process like in 
the case r = 2 treated in the beginning of the proof. Then the function / will have to satisfy 



/(g) > d 2 / m (min{g, r}) + ^ dA(min{g, r}), 



i=3 



with some dj > for all i > 2. Eventually, e must be small enough such that f(q)e < 1. But / 
is increasing with q and fmie f fixed, there exists a rank k = k{e) where the above requirement is not 
satisfied. Eventually, Th. ll4.ll gives |Pi| = e and this concludes the proof. □ 



stable 



6 Bifurcation point and stability 

In this section, the regime of the system around the bifurcation point is investigated. The direction of 
the bifurcation, the amplitude of the first Fourier component and the influence of 7 on the oscillation 
frequency are given. The results are similar to the ones arising in the finite dimensional Hopf bifurcation. 
As for the Hopf theorem, an additional hypothesis on the system is required, i.e. the transversality 
condition stated in H4. 

Without loss of generality, as explained in the previous part, the index e of H2 is assumed to be e = 1 
in the following. Suppose H2 holds, then for 7 6 f2, /^(t), the eigenvalue of J 7 associated to u>2, is not 
null. Assume H3 and let us introduce the notation: 

£1(7) = [foil)]- 1 92(2, 1) (, 92 (1, 2) + 0,(1, -1)) - (53(1, 0, 1) + 0,(1, 0, -1) + g a (l, 2, 1)) . (50) [drfD 

Let us denote by lZe[z] the real part and Im[z\ the imaginary part of the complex number z. In this 
section, the following additional hypothesis is assumed: 

H4 There exists an open set T with 70 S T such that the complex-valued functions D\ and (3\ are 
respectively continuously differentiable and twice continuously differentiable with respect to 7 on 
r. Moreover, -Di(7o) 7^ and 

-Di(7o) «7 
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ThPl | Theorem 6.1. Suppose HI, H2, H3, H4 hold. Then 
i) The quantity 

1 d/3 



-Di(7o) ^7 



-(70), (51) def alpha 



is a real number, different from zero. If its sign is positive (resp. negative) the bifurcation is direct 
(resp. inverse), i.e. the oscillating solution described in HI exists for 7 > 70 (resp. 7 < 70J. 

ii) for all 7 S T such that I7 — 7o| is small enough, the modulus of the first Fourier component of the 
oscillating solution is given by: 



\Pi\ = VW\VW^\+0(\7~Jo\ 3/2 )- (52) 

JTM I I firnjnflSI 

Proof. For e small enough, from Th. lb. II and lemma IB. 41 one can write: 

S(g 2 , 2, 1) = (g 2 {l, 2) + g 2 (l, -1)) P^P 2 + 0(e 4 ), 

S(g 3 , 3, 1) = ( 53 (1,0, 1) + . 93 (1, 0, -1) + g 3 (l, 2, 1)) \P 1 \ 2 P 1 + 0(e 4 ), 

Vn > 4, |S(s„,n, 1)| < cie 4 , with ci > 0. 

One has also: 

5( ff2 ,2,2)= 52 (2,l)P_ 1 P 1 +0(e 4 ), 
Vn>3, \S(g n ,n,2)\ < c 2 s\ 
Vn>2,g>3, |5(ff„, n, q)\ < c 3 e 4 . 

Then it yields for P 2 : 

00 

J^2 = /5 2 (7)^2 = 51 2 ) = 32(2, 1)P-iPl + 0(e 4 ) 

i=2 

p 2 = [ftwr '52(2, i)p-iPi + o( £ 4 ), 

For Pi: 

OO 

J 7 p i = ^2 s (9t,h 1) 

i=2 

= (g 2 (l, 2) + 52 (1, -1)) P_iP 2 + («? 3 (1, 0, 1) + 53(1,0, -1) + 53(1, 2, 1)) |Pi| 2 Pi + 0(e 4 ), 



which yields 

Pi{/3 1 ( 7 )-Pi(7)| J Pi| 2 }+O(e 4 ) = 0, 

where D\ is the differentiable function defined in (|bu|). For the oscillating solution, where Pi 7^ 0, one 
can write 

|P 1 | 2 Pi( 7 )=/3 1 ( 7 ) + 0(e 4 ), 
From the hypotheses of H4 and H2, one has around 70 the following Taylor expansions: 
Px( 7 ) = Pi (70) + 0{Yi - 7o |), ft (7) = #( 7 o)(7 - 7o) + 0(7 - 7o) 2 - 

and 

|Pi f = ~ 7o) + 0(7 ~ 7o) 2 + 0{e% (53) 

£>i(7o) 

JThhi f 

For e | grn all enough, one has |Pi| = e from Th 114.11 and since the numerator is not nulljxorn H4, equa- 
tion [p'6\) shows that (7 — 70) is proportional to e 1 . Moreover, the left hand side of (|ba|) is a positive 
quantity, implying: 

7ee[a](7-7o)>0 (54) 
Zm[a\ = 

Hence the sign of 72.e[a] determines the direction of the bifurcation (direct or inverse). □ 



eqPl 
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The last result of this section concerns the influence of 7 on the oscillation frequency, u)\, around 
the threshold. It requires two further hypotheses, one on the linear part of the system and one on the 
coefficients in front of the quadratic and cubic terms, contained in the function D\ . The proof is in fact 
an application of the implicit function theorem. For the next result, (3\, D\ shall be written as functions 
of two variables: /3(wi,7) = D(ivi,-f) = Di(j). 



corolomega Corollary 6.2. Under the hypotheses of Th. lp.L\) . suppose the function Tm\D(3] is twice continuously 
differentiable with respect to its variable lo\ and 7 in an open set around the point (oji,7q), and 



^0, 



("1,70) 



doj-f' 



—ImlDfc] 



?0. 



Suppose further that [3 is differentiable with respect to both variables u>\ and 7 at (0^1,70). Then the 
expression of the angular frequency uj\ for 7 in a neighbourhood of 70 is given by: 



(55) 



^1(7) = ^1(70) +^i(7o) • (7 - 7o) + 0(7 - 7o) 2 , 
where (^1(70) is the solution o/Xm[/3(ti;i(7o), 70)] = and 

^i(7o) = 
J amegaT 

Proof. The Taylor expansion in Eq. (|bij|) is a consequence of the implicit function th qp^gr n applied to the 
twice differentiable function IrrdDB] around 70. The value of cji(7o) is given by Th. fe.il The hypothesis 
H2 2) and the point i) of Th.inT imply that 



Im 




(wi(7o),7o) 


Tm 


OUJ 


(^i(7o),7o) 



Xm \ D l X ](7o) =Im £»i(7o)^-(7o) 
dj 07 



= 0. 



(56) 



The implicit function theorem and the hypotheses of differentiability allow to write 



Im 



D 



dpi 
f/7 



(7o) 



djj 
f/7 



(70) • Im 



D 



(wi(7o),7o) + Tm 



D 



97 



(wi(7o),7o)- 



Putting this latter expression in Eq. 



dIO 

H|) concludes the proof. 



□ 



omegaT 



condIO 



woodwind 



7 Example of the woodwind musical instruments 

A reed woodwind musical instrument, such as a clarinet or a saxophone, is an example of a system which 
can be described by an hyperbolic PDE and a localized polynomial nonlincarity. In acoustics, the major 
challenge is to understand the influence on the produced sound of physical parameters such as the pressure 
in the mouth of the musician, the stiffness of the reed, the shape of the mouthpiece and of the resonator. 
Assuming it can be done, a first order calculation by linearizing the system gives the bifurcation point and 
the frequency of oscillation at this point. But the investigation of these systems requires to go further 
and to obtain the shape of the Fourier harmonics with respect to the bifurcation parameter. That is 
why the frequency domain approach is perfectly appropriate here. For small oscillations, one can obtain 
analytical formulae giving information on the solution of the system, but, as said in the introduction this 
requires until now a hypothesis on the Fourier series of the solution. Let po £ L 2 (T) be the acoustic 
pressure inside the mouthpiece of the wind i : 9^'TOi cn ^ K r /r a dS r (g^ uc ly ana ^ ^ us suppose it is oscillating. 
The hypothesis made by acousticians, see e.g. \27 , [TB"~ (T7~s that the g-th harmonic of 
obeys: 



: of this pressure 



\P q \<\Pl\ M , q^O. 



(57) 



eqWorman 



This hypothesis is now going to be proved in the framework of the mathematical results obtained in the 
previous sections. _ 

The equations describing this system are detailed in e.g. |fj, 0. Although the instrument possesses 
a three dimensional shape, the model is one dimensional and X = [0,£] where t is the length of the 
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instrument. The dimcnsionless variables associated to the system are the acoustic pressure, p, and the 
volume flow, u, inside the instrument, both belonging to X, and h <E L 2 {T) the reed tip opening localized 
x = 0. The pressure inside the mouth of the instrumentalist 7 is the bifurcation parameter. Let us call 
Po and uq the pressure and volume flow at x = 0. The acoustic pressure depends on the volume flow 
entering the mouthpiece through the relation: 

Ap = 5{x)Bu, (58) 

where B = d/dt and A is the differential operator associated to the wave equation inside the instrument, 
with losses at the boundaries and along the resonator (visco-thermal losses) . By using the Green functions 
formalism, it can be written as 

p(t,x) = g(t-T,x-y)Bu(T,y)S{y)dydT = g(t~T,x)Bu(T,0)dT, (59) 

Jtxx Jt 

where g is the Green function solution of Ap = 6(x)5(t). An example of this Green function for a model 
of clarinet is given in n?, Chapter II. 5, section 5.2: the model considered is the wave equation with 
visco-thermal losses on X, Neumann boundary condition at one end and Dirichlet boundary condition 
at the other end. This gives for the time Fourier transform 7j of g: 

s-^ cos(k n x)cos(k n x s ) 

g{x,u;x s ) = c > 60 

^ w « + iuu n q n - UJ Z 

where c is a constant, k n = (2n — 1)tt/(2£) and w„ is proportional to k n . The term q n is proportional to 
\/fcn in order to model visco-thermal losses. This choice of q n is motivated by physical arguments but 
leadsjx^a^fractional derivative in the expression of A. This is partly why the present study start with 
Eq. (jlj, where Z is much more convenient to define than A. Let U = J-uo, the above equation expressed 
in the Fourier space, and at x — 0, yields: 

P = ZU, (61) 

where Z is called a transfer function. This is the equation (jH]) m the case of wind instruments, in the 
frequency domain. Indeed, U is expressed in terms of P through a nonlinear relation in the following. 
Notice that the quantity Z is called by acousticians the input impedance of the bore and can be measured 
thanks to physical experiments. It is bounded, provided losses are taken into account in the model. See 10, 
chapter II. 7, for examples of impedances of various wind instruments. The Fourier basis is the cigenbasis 
of this operator, thus in the Fourier space, Z is a multiplication operator. For all q € Z, 

P q = Z(uj q )U q , (62) [PZU2" 

where the complex-valued function Z(-) is never null except maybe at luq = 0. 

The nonlinear relationship is localized at the tip of the mouthpiece and relates uq, po and h. Under 
some physical hypotheses, the Bernoulli law at the entry of the instrument gives the nonlinear relationship 
at time fel: 

ug(f) = ( 2 (1 - + h(t) f ( 7 - p (t)) , (63) I eqnlt2 

where £ is a positive constant, < 7 < 1. At rest, when there is no action from the instrumentalist, 
7 = 0, and the equilibrium position of the instrument yields uq = 0, po = and h = 0. It is supposed that 
for all t, po(t) and h{t) are small enough such that 7 — po(t) > and 1 — 7 + h(t) > (small oscillations). 

Let H denote the Fourier transform of h. The reed opening is related to p$ via a linear relation 
expressed in the frequency domain: 

H q = D(q)P q , for all q G Z. (64) | EDP 

The complex- valued function D is bounded. Let us introduce Qq the orthogonal projection associated to 
the constant part of the Fourier decomposition. From ip'Z]) . one has 

U q = Y(q)(l-Q )P q + U Q Q 0l forallqeZ. (65) [pzu 
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In the physical system under study, the function Y(l — Qq) : Z — * C is bounded. Nevertheless, the value 
Y(0) is not defined as it is assumed that Pq = and Uq > 0. It gives a good opportunity to show how 
to cope with this problem and satisfy the hypothesis of boundedness of the operators despite it: in this 
case the^vstem wlUJbe 2 written as a set of two equations, with a specific one for Uq, see below how to 
38 b|) . Behind the function Y is hidden the hyperbolicity given by the wave equation 



obtain (|K8a|) and 

inside the resonator and the contribution of B. It reveals itself when looking at the limit of |5^(?)| when 
q tends to infinity: it tends to a strictly positive constant. 

Let us assume the pressure po satisfies hypothesis HI. Let us recall that since D is bounded, one has 
for example, for all geZ, 

hp(q) = (H * P)(q) = ((DP) * P) (?) = Dg- n P q - n P n , 

nGZ 

naiTitimtf _ _ IRGKSV 

as shown by lemma IB. II It. like in Ijl Yj. the following notations are introduced: 



eqNLP2T 



B 



q.n 



C a 



C 2 7 (l - 7 ) 2 , A q = 2C 2 7 (1 - l)D{q) - C 2 (l - 7) 2 
= C 2 lD(q - n)D(n) - 2( 2 (1 - j)D(q - n), 
„ = -( 2 D(q - n)D(n - m), 



and for n 7^ 0, n 7^ q, 



n q ,n = B q ,n ~ Y(q - n)Y(n), 



then the relation 



lt2 

TTbgether with 



RGKSV 

gives (see I|1Y] tor more details): 



(2U Y{q) - A q )P q = J2 K q ,nP q -nPn 



EE 



Cg,n,mPg- 



P P 



n^Z m£Z 



P P P 



for q f 
for q = 0. 



(66) 
(67) 

(68a) 
(68b) 



Qo) and D are bounded and satisfy H3. As said previously, 

'"--" reg 

one has 



Lemm a lb . 1 1 1 ustifics this expression as Y(l 

Eq. (|b8b|) appeaJs jjeca use Y is not defined a tO^fa order to obtain an expression similar to |T2 
to replace Uo in (|b^a[) by its value given in (|bab|) . For this purpose it is assumed that Uq has a positive 
value, allowing to express it as a square root. It is assumed that the system encounters a bifurcation 

pressure, this square root can be 



for 70 > which implies uoo > 0. Thus, for a small enough acqu: 
expanded as a Taylorserie|. Replacing Uo by this expression i r^jpSal 
presented in lemma B.li As a consequence, theorems fe.il and fe.il can be used 



complex to write in the canonical form given in 

J 7 (?) 




leads to a nonlinearity of the form 
The above relations are 



"however, for instance one has for the linear part: 



U oY{q) - A q ). 



eqNLP2 



eqNLP20 



(69) eqJcar 



This linear part satisfies H2, provided (2^JuooY(q) — A q ) = for only two values q — ±s. Notice that 
the stationary solutions Uq = ±^/uoo, where for all n E Z, P n = 0, exist for all 7 G (0, 1). Consequently, 
from (|Kij)|) . trie bifurcation parameter and the frequency of oscillation at the bifurcation will be given, 
according to Th. [OJ by the characteristic equation: 

- A s = 0. 



The main objective is now attained: hypgthescs^ H2 and H3 hold and if one supposes a solution of 
the type described in HI, the hypothesis (|py|) can be justified in the framework of the present study by 
applying Th. j^J ^Thc results presented in the publications cited in this section, which required at first the 
assumption (|b 7 1) . can be obtained by assuming the weaker assumption of twice continuous differentiabUifcy 
of the solution. Let us notice a slight difference between, the result of theorem lb. II and the relation (lb 7 1) . 
The fact that there exists a value k in Th. lb. II where (|b?|) stops to be valid does not change the results of 
the previously cited publications. Indeed, these results can be retrieved by knowing estimates on the first 
harmonics and involve only a limited number of harmonics in calculus. Moreover, k can be taken as large 
as needed by limiting the range of amplitude for which the results are valid. Eventually, the hypotheses 
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required in section rr are satisfied, see l[TV|. and the quantity -Di(7o) defined in this latter section can be 
deduced from equations fooji . This leads to 

\Pi\ = vHvt^to + o(h - 70 | 3/2 ), 

where 

= £M_ 

" D 1 (70)' 

IRGKSV 

confirming the results of [T7fr 
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